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ABSTRACT 

Infrared vibration-rotation lines can be valuable probes of interstellar and circumstellar 
molecules, especially symmetric molecules, which have no pure rotational transitions. But most 
such observations have been interpreted with an isothermal absorbing slab model, which leaves 
out important radiative transfer and molecular excitation effects. A more realistic non-LTE and 
non-isothermal radiative transfer model has been constructed. The results of this model are in 
much better agreement with the observations, including cases where lines in one branch of a 
vibration-rotation band are in absorption and another in emission. In general, conclusions based 
on the isothermal absorbing slab model can be very misleading, but the assumption of LTE may 
not lead to such large errors, particularly if the radiation field temperature is close to the gas 
temperature. 

Subject headings: Line:formation - Radiative transfer - ISM:molecules 



1. The Promise 

Infrared spectroscopy of molecular vibration- 
rotation lines holds the promise of providing a 
valuable probe of chemical and physical conditions 
in interstellar and circumstellar gas. Vibration- 
rotation spectra complement pure rotational spec- 
tra in several important ways. First, molecules 
without permanent dipole moments, such as CH4, 
C2H2, and CO2, can be observed when vibration 
breaks the symmetry of the molecule. Second, be- 
cause vibration-rotation lines are typically seen in 
absorption toward embedded young stars they se- 
lectively probe gas in star-forming regions. And 
third, because vibration-rotation lines of different 
rotational states lie close together in wavelength 
they can be observed with one instrument and 
telescope, and often simultaneously. 

It seems that the interpretation of absorption 
lines should be straightforward. The optical depth 
of an absorption line depends primarily on the col- 
umn density in the lower state, that is Nj(v = 0), 
so the calculation of Nj from the spectrum is sim- 
ply a matter dividing the equivalent width of a 



line by the line-strength factor obtained from lab- 
oratory spectroscopy. If many lines are observed, 
it is not even necessary to assume a thermal dis- 
tribution of lower state populations, as the pop- 
ulation of each rotational state can be measured. 
Of course, for optically thick lines curve-of-growth 
or other corrections for saturation may be neces- 
sary. This introduces uncertainties due to the gen- 
erally unknown lineshape, but the problem can be 
avoided if necessary by observing optically thin 
lines of rare isotopes. This interpretation proce- 
dure is referred to as the isothermal absorbing slab 
model. Variations on it have been used by as- 
tronomers to interpret almost all infrared molecu- 
lar absorption spectra. 

2. The Problems 

There are several problems with the isothermal 
absorbing slab model: the dust which emits the 
'background' continuum radiation is often mixed 
with the absorbing gas, the molecular gas emits as 
well as absorbs radiation, and there is a tempera- 
ture gradient in the gas and dust around embed- 
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ded sources. 

As a first approximation, tlie effect of mixed gas 
and dust can be taken into account by consider- 
ing the derived molecular column densities to be 
the columns to a depth in the source where the 
dust optical depth, r, reaches one. But since the 
dust opacity is wavelength dependent, this depth 
generally differs for different molecules or differ- 
ent vibration-rotation bands, and may even dif- 
fer substantially for lines within a band. In addi- 
tion, the emergent continuum radiation at a given 
wavelength is not necessarily emitted from near 
T = 1. If there is a temperature gradient in 
the dust, with the dust in the source being hot- 
ter at greater depth, emission by dust at r > 1 
may dominate the outgoing radiation, especially 
at wavelengths shortward of the peak of the Planck 
function. Without a radiative transfer calculation 
that takes these effects into account, the column of 
gas through which absorption line observations are 
made may be substantially underestimated, and 
even derived molecular abundance ratios may be 
in error. 

The effect of emission in the molecular lines is 
even more difficult to determine and account for. 
If a source is spatially resolved and the molecular 
gas is in LTE at a known temperature, the equa- 
tion of radiative transfer along the line of sight is 
easily solved. The result is that molecular lines 
saturate at an intensity equal to the Planck func- 
tion at the line wavelength. But often the contin- 
uum source is unresolved and the absorbing and 
emitting gas may have different angular extents, 
which may not fall entirely within the observing 
beam. Even more problematical is the fact that 
collisional deexcitation cross sections are small 



([Gonzalez- Alfonso et al.ll2002r ) and radiative rates 
are relatively large for vibration-rotation transi- 
tions, making the critical densities for vibrational 
thermalization large, typically > 10^^ cm~^. Con- 
sequently, vibrational LTE at the gas kinetic tem- 
perature could be a very poor approximation. It 
may be a better approximation to assume that all 
radiative excitations are followed by radiative de- 
cay. In a spherically symmetric situation in which 
the molecules lie in a shell separated from the con- 
tinuum source, but within the observing beam, 
there may be no net absorption even if the lines 
are optically thick, since an equal number of pho- 
tons are emitted into the beam as are absorbed out 



of the line of sight. The actual net absorption or 
emission depends on whether more or less molec- 
ular gas lies along the observing line of sight than 
along a typical direction and whether reemitted 
photons are absorbed by dust before escaping from 
the source region. In the case of sources with sys- 
tematic motions, notably expanding shells around 
evolved stars, the emission and absorption may 
be separated spectrally by the Doppler effect. In 
those cases high spectral resolution observations 
would avoid this problem, but in many other cases 
the systematic motions are smaller than turbulent 
linewidths or are unresolved, making it impossible 
to know how much the emission and absorption 
cancel. 

In many cases the radiative transfer effects we 
have been discussing change the strengths of lines, 
so lead to errors in derived abundances, but are 
not obvious from inspection of the spectra. But 
in a few cases there are recognizable symptoms 
of these effects. The importance of emission in 
vibration-rotation lines is clearest in the spectra 
of expanding circumstellar shells, with IRC-|~10216 
being a particularly good example. The bending- 
mode Q branches of HCN and C2H2 near 14 /xm 
are prominent in t he R=2000 ISO SWS spec- 
trum of IRC+10216 (jCernicharo et al.lll999t ). with 
depths of ~ 20% of the continuum, but the P and 
R branch regions of the spectrum show little evi- 
dence of line absorption or depression of the con- 
tinuum due to lines. In contrast, the R=80,000 



TEXES observations of iFonfrfa et al.l (|2008l ) show 
the R branch lines of HCN and C2H2 to have 
typical depths of ~ 50%. The reason the lines 
are so much weaker in the ISO SWS spectrum is 
that the lines have P-Cygni profiles, with nearly 
equal emission and absorption components that 
are blended together in the R=2000 spectrum. 
Apparently most vibrational excitations caused by 
absorption of photons in these bands are followed 
by emission, and most emitted photons escape 
from the shell, perhaps after multiple absorption 
and emission events. 

The 13.7 /im 1^5 band of C2II2 has also been ob- 
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Ba rentine fc Lacvl ()2012n also observed absorption 
in the 7.5 /im i'4+i'5 band toward OMCl IRc2 and 
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Fig. 1. — C2H2 vibrational energy levels and tran- 
sitions of interest. The 1^4 and 1^5 modes are the 
anti-symmetric and symmetric bending modes, re- 
spectively. Allowed radiative transitions are indi- 
cated with solid arrows; the Vi transition that only 
occurs coUisionally is indicated with a dashed ar- 
row. Allowed radiative transitions between these 
levels must change the quantum number V5; V4 
can change only if V5 also changes. The 85 /im 
1^5 — V4 transitions have not been observed. Ro- 
tational levels and level splitting are not shown. 
The vertical scale is in cm~^. 



NGC 7538 IRS 9. An energy level diagram of C2H2 
with the relevant vibrational levels is shown in 
Figure 1. Both bands involve absorption from the 
ground vibrational level, but the derived C2H2 
column density from the 1^4-1-1^5 lines is greater 
than that derived from the lines by a factor of 
3-10 toward IRc2 and 20 toward IRS 9. Both of 
the effects described above may contribute to this 
discrepancy. The dust opacity at 13.7 /im is about 
twice as large as that at 7.5 /im, because of the 
wings of the 9.7 and 18 fj,m silicate features. In 
addition, 7.5 /zm is farther off of the peak of the 
dust thermal emission, also causing the continuum 
to be formed farther into the sources (where the 
temperature is higher) at 7.5 /xm. The amount 
of reemission in spectral lines may also differ be- 
tween the 1/5 and the 1^4+1^5 bands, since the 
state can only decay to the ground via the 1/5 
band, whereas the U4+U5 state can decay either to 
the ground through emission of a 1^4-1-1^5 photon or 
to the 1^4 level by emission of a i'4+iy5-i>4 photon. 

Ammonia (NH3) als o shows evidence of ra- 
diative transfer effects. iBarentine fc Lacy ( 2012 ) 
observed lines in two branches of the NH3 1/2 
band toward NGC 7538 IRS 9. Although other 
molecules seen toward IRS 9, including C2H2 and 
HCN, are seen in absorption and the aQ branch 
of NII3 shows P-Cygni lines or weak absorption, 
aP and sP-branch NH^, lines are seen in emission. 



Barenti ne fc Lacvl (j2012n briefly discuss a radia- 
tive transfer model to explain these observations. 
It assumes that the NII3 molecules are exposed to 
optically thin silicate emission, which is stronger 
on the NH3 R and Q branches, so causes more 
upward transitions in those branches than in the 
P branch. Since downward transitions are about 
equally likely in the three branches, emission dom- 
inates over absorption in the P branch, and ab- 
sorption dominates over emission in the Q and R 
branches. 

3. A Radiative Transfer Model 

To reproduce and understand the radiative 
transfer and radiative excitation effects in molec- 
ular vibration-rotation spectra, I constructed a 
computer model with includes collisional and ra- 
diative excitation and deexcitation of molecules 
and calculates the outgoing spectrum. 

The model assumes a parameterized distribu- 
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tion and composition of gas and dust in a shell 
around a blackbody luminosity source. The input 
parameters can be adjusted through a grid search 
procedure to fit the continuum and line spectra of 
individual sources. For the calculations in this pa- 
per we assume a spherically symmetric expanding 
shell. The model first calculates the dust temper- 
ature and continuum radiation field by alternately 
solving the equation of radiative transfer and the 
equation of thermal equilibrium for the dust. The 
radiative transfer calculation involves following a 
grid of rays through the shell. Typically, a 64 x 64 
ray grid is used at the outer edge of the model. Af- 
ter following the rays half way in to the star, the 
inner 32x32 rays are split into 64x64 rays to in- 
crease the resolution. This is done up to 11 times, 
and then the rays are rejoined going out from the 
star. Rays that overlap the star are replaced by 
an appropriate fraction of the stellar spectrum be- 
fore continuing out through the shell. As the rays 
are followed through the shell, the ray spectrum 
is added to a mean intensity spectrum stored for 
each grid cell, which is used to calculate the dust 
heating. The dust temperatures are calculated for 
cells on a spherical coordinate grid. For each grain 
type and size the equilibrium temperature is calcu- 
lated that balances radiative heating and cooling. 
The dust source function is then recalculated and 
the radiative transfer calculation repeated. Typ- 
ically 16 iterations of the radiative transfer and 
thermal equilibrium calculations are sufficient for 
even very optically thick dust shell models to con- 
verge. The gas temperature is not calculated sep- 
arately, but is assumed to equal the dust temper- 
ature. For the models described in this paper, the 
dust shell and stellar parameters were chosen to 
roughly reproduce the observed SED of NGC 7538 
IRS 9. 

After solving for the shell thermal structure, 
the program solves for the molecular level popu- 
lations and the outgoing line spectrum. This is 
done by alternately solving the equation of ra- 
diative transfer in molecular lines and the equa- 
tion of statistical equilibrium involving collisional 
and radiative rates in and out of vibrationally 
excited energy levels. Line strengths are calcu- 
lated from band strengths an d Honl-London fac- 
tors, as in lEvans et al.l (jl99l[ ). and are consistent 
with those in the HITRAN database. The calcu- 
lation is simplified for linear and symmetric top 



molecules by the selection rules that cause each 
vibrationally excited rotational level to be cou- 
pled radiatively to only two or three rotational 
levels in the ground vibrational state. It is fur- 
ther simplified by assuming that the rotational 
levels in the ground vibrational state are in LTE 
at the gas temperature. This and other assump- 
tions are discussed in §5. We also assume that 
collisional transitions between vibrational states 
do not change the rotational state. This is not 
actually valid, but since collisions must drive pop- 
ulations toward LTE, and the ground vibrational 
state is assumed to be in LTE, on average col- 
lisions should not change the rotational excita- 
tion of molecules. In addition, vibration-changing 
collisional transitions are relatively unimportant 
compared to radiative transitions at the densi- 
ties found in the model. With these simplifica- 
tions, the excited vibration-rotation level popu- 
lations and the vibration-rotation lines between 
them and the ground vibrational state can be cal- 
culated for one excited vibration-rotation level at 
a time, neglecting any coupling between rotational 
levels of the excited vibrational state. 

As with the dust calculation, the code follows 
a grid of rays through the shell, in this case calcu- 
lating the spectrum for wavelengths close to the 
relevant vibration-rotation lines. At each step 
through the shell the molecular line profile is 
Doppler shifted by the component of any assumed 
systematic motion along the ray direction. As 
the rays propagate through the shell the Doppler- 
shifted spectrum is added to a sum for the rays 
passing through each spherical coordinate cell, so 
that the mean intensity, averaged over each spec- 
tral line, can be calculated for each cell. With the 
mean intensities, the radiative excitation and de- 
excitation rates are calculated, and with those and 
the collisional rates, the vibration-rotation level 
populations are calculated. The radiative trans- 
fer and level population calculations are then re- 
peated until they converge. However, the fact 
that the molecular lines are typically more op- 
tically thick than the dust continuum results in 
slower convergence than for the calculation of the 
thermal structure. This problem was overcome 
by comparing the change in level populations in 
each iteration with that in the previous itera- 
tion, and amplifying the correction if iterations 
are well correlated. With the accelerated conver- 
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gence scheme, the level populations almost always 
converge within 16 iterations. 



4. Results of Modeling 



Models were run with parameters meant to 
roughly reproduce the IRS 9 observations. The 
central source was taken to be a 12,000 K, 10^ L0 
blackbody. The gas and dust shell had a con- 
stant H2+He density of 2 x 10''cm~^ from 4 to 
400 AU, then fell with an density profile out 
to 1.6 X 10** AU. Its mass was 4.2 M©, and its ra- 
dial column density was 1.8 x 10^^ cm~^. The cal- 
culated dust temperature was 1800 K at the inner 
edge of the shell, falling to 160 K at 400 AU and 
40 K at the outer edge of the shell. A constant ex- 
pansion speed of 2kms~^ and a turbulent veloc- 
ity dispersion of 3kms~^ were assumed. Molec- 
ular line profiles were calculated for the C2H2 1/5 
and 1^4+1^5 bands, with C2H2/II2 = 1 x 10"'^, the 
NH3 1^2 band, with NH3/H2 = 4 x 10"^, and the 
CO v=l-0 band, with i3CO/H2 = 1 x 10"^ and 
C*^0/Il2 = 2 X 10"'^. These parameters are not 
meant to be fitted values, but they produced line 
depths similar to those observed. 

Several of the calculated NH3 line profiles are 
shown in Figure 2. In each panel, the P, Q, and 
R-branch lines between a given J,K state in the ex- 
cited vibrational level and the ground vibrational 
level are shown. Note that these are normally re- 
ferred to as the P(J„-|-1,K), Q(J„,K), and R(Jm- 
1,K) lines, being labeled by their rotational levels 
in the lower vibrational state. For this 'umbrella' 
vibrational mode, Q-branch lines have the largest 
optical depths for K=J, and P and R-branch lines 
are strongest for K=0. The nuclear statistical 
weights favor lines with K=3n by a factor of 2. 
In all cases, the model P-branch lines are pre- 
dominantly in emission, the Q-branch lines show 
weak absorption, and the R-branch lines are more 
strongly in absorption. This is in agreement with 
the observed P and Q-branch lines toward IRS 9. 
The R-branch lines have not been observed, but 
are predicted to be seen in absorption, with depths 
of 1-2%. 

Several calculated C2II2 ly^ lines are shown in 
Figure 3, and va+vq lines are shown in Figure 4. 
For C2H2, rotational states in the ground vibra- 
tional level with Jj odd (ortho states) have nuclear 
statistical weights of 3; those with J; even (para 




-40-20 20 40 -40-20 20 40 -40-20 20 40 



Fig. 2. — Model NH3 V2 umbrella-mode line pro- 
files. In each panel, the P, Q, and R-branch lines 
between a given J,K state in the excited vibra- 
tional level and the ground are shown (top to bot- 
tom). Spectra are normalized, with the Q and 
R-branch lines offset downward. Q-branch lines 
are strongest for K„ — J„ and are absent for K„ 
= 0. R-branch lines are absent for = Z^. 
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states) have nuclear statistical weights of 1. Line- 
strength factors for P, Q, and R-branch lines are 
proportional to J„, 2J„+1, and J^+l, respectively. 
Combining these factors for the lines shown, all of 
which have J„ even, so J; odd for the P and R- 
branch lines, the P and R-branch lines have sim- 
ilar line-strength factors, whereas Q-branch lines 
are weaker on average by a factor of 2/3. In addi- 
tion, for absorption lines the line-strength factors 
must be multiplied by the Boltzmann population 
factors of the lower rotational states. 

In agreement with the IRS 9 observations, the 
U5 R-branch lines have depths only about twice 
those of the Vi+u^ lines, even though the band 
strength is nearly 10 times greater. The P-branch 
line depths are similar in the two bands. A new 
prediction is that the 1/5 P-branch lines should be 
substantially weaker than the R-branch lines. At 
higher abundances the P-branch lines go into emis- 
sion and the Q-branch lines show P-Cygni pro- 
files. The P-branch lines are unobservable from 
the ground, due to telluric CO2 absorption. The 
Q-branch lines are difficult to observe, due to tel- 
luric absorption and blending of the closely spaced 
lines, but do show some evidence of P-Cygni emis- 
sion. The explanation of the strengths and profiles 
of these various lines is discussed below. 

Calculated CO lines are shown in Figure 5. CO 
has no Q branch; the third line shown in each 
panel is for an average of the P and R-branch line 
strengths, but with emission by molecules in the 
v=l state turned off in the program, to show what 
would be observed if the pure-absorption model 
used in the past were valid. Both P and R-branch 
lines of CO have P-Cygni profiles, with the rel- 
ative strength of tlie emission component being 
greater in the P-branch lines. The observed lines 
of CO from IRS 9 also have P-Cygni profiles, but 
are considerably broader than the lines of other 
molecules. This is probably a result of a high ve- 
locity outflow with a different chemistry from the 
gas that dominates the spectra of other molecules. 
The high optical depths of the ^^CO and ^^CO 
lines may also contribute to the prominence of line 
wings that are not apparent in spectra of otlier 
molecules. C^^O and C^'^O lines are blended with 
^^CO and ^^CO lines in the existing observations; 
further observations would be desirable. 
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Fig. 3. — Model C2II2 symmetric bending mode 
line profiles. In each panel, the P, Q, and R-branch 
lines between a given J„ level in the excited vi- 
brational level and the ground are shown (top to 
bottom). Q and R-branch lines are offset down- 
ward. Rotational levels of the excited vibrational 
level are split into ortho and para states, with the 
P and R-branch lines going to the ortho (g„ = 3) 
state for J„ even, and the Q-branch line going to 
the para (g„ = 1) state. 
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Fig. 4. — Model C2H2 1^4+1^5 band line profiles. 
Only P-branch (top) and R-branch (bottom) lines 
occur for this band; the middle spectrum is for 
the 1^4 + f5 — 1^4 Q-branch line, which accounts for 
much of the radiative deexcitation of the 1^4+1^5 
state. Since it is scaled by the continuum flux, the 
I'A + 1^5 — line appears weaker than the v^+vq 
lines. 



Fig. 5.— Model ^^CO v=l-0 line profiles. Only 

P (solid) and R-branch (dotted) lines occur for 
diatomic molecules like CO. The dashed line 
in each panel is for an average of the P and R- 
branch line strengths, but leaving out emission in 
the line. 
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5. Explanation and Discussion 

The model spectra are generally in at least qual- 
itative agreement with the observations of IRS 9. 
In particular, they reproduce the emission seen 
in NH3 P-branch lines and the similar depths of 
C2II2 1^5 and Vi+vs lines. However, they make 
additional predictions not antic ipated by the sim- 
ple ra diative transfer models of iBarentine &: Lacv 
(|2012 h. Notably, they predict that C2II2 should 
show differences between its P, Q, and R-branch 
lines like those seen in NH3 . Somewhat less promi- 
nent effec ts are predicted for CO. T he explanation 
given by IBarentine &: Lacvl ( 2012 ) for NH3 was 
based on the brighter continuum radiation at the 
Q and R-branch wavelengths. But this should not 
be the case for C2H2, for which the dust opacity is 
similar at the wavelengths of its different branches. 

To help understand the radiative transfer ef- 
fects responsible for the line profiles of different 
lines, several physical effects were changed in the 
model. First, a model was run assuming LTE 
populations of the excited vibration-rotation lev- 
els. Perhaps surprisingly, the model results did not 
change greatly; P-branch lines still favored emis- 
sion relative to R-branch lines. Apparently the 
non-LTE populations were not very different from 
LTE populations. The explanation for this may 
be that the continuum radiation field seen by the 
molecules is not very different from a blackbody 
field at the dust temperature, and the gas temper- 
ature was assumed to be equal to the dust tem- 
perature. As a result, populations controlled by 
radiation were not very different from populations 
controlled by collisions. This would not be the 
case if the molecular gas were more distant from 
the warm dust responsible for the continuum, so 
the conclusion that non-LTE effects were relatively 
small for the IRS 9 model may not be valid in gen- 
eral. A second modification to the physics of the 
model was to set the rotational constant (Bg) of 
the molecules to zero, so that the P, Q, and R- 
branch lines all fell at the same wavelength, and 
their lower state energies were all the same. The 
differences between P, Q, and R-branch lines then 
disappeared. Finally, a model was run with the 
gas and dust temperature held constant through- 
out the shell. Again, the differences between P, Q, 
and R-branch lines largely disappeared. The con- 
clusion from these models with modified physics 



was that the difference between the lower state en- 
ergies, and the associated Boltzmann factors, was 
the primary cause of the different line profiles. 

With the results of the modified models, the 
explanation for the differences between the P, Q, 
and R-branch lines can be understood. We will 
consider P and R-branch lines, which have simi- 
lar behaviors for all molecules. (The strengths of 
Q-branch lines depend on the molecular symme- 
try and vibrational mode.) The P and R-branch 
lines from a given rotational state of the upper vi- 
brational level have nearly equal Einstein A coef- 
ficients (depending weakly on J„), so nearly equal 
branching ratios on emission. However, the ab- 
sorption in these lines depends on the radiation 
field and the populations of the rotational states of 
the lower vibrational level. Especially for high val- 
ues of J, these populations can differ substantially, 
because of the different energy levels of J; — Ju + i 
and Ji — Ju — I for the P and R-branch lines, re- 
spectively. As a result, the two lines have simi- 
lar emission strengths, but the R-branch line has 
greater absorption. If most upward transitions are 
followed by a downward radiative transition, emis- 
sion and absorption nearly cancel, but with emis- 
sion dominating in the P-branch line and absorp- 
tion dominating in the R-branch line. In the case 
of NH3 the greater radiation field at the wave- 
lengths of the R-branch lines enhances the effect, 
but the effect is present even for C2H2 and CO in 
models with a temperature gradient through the 
shell. 

The disappearance of the difference between P 
and R-branch lines in the isothermal model can 
be explained by the fact that the continuum radi- 
ation field at the R-branch wavelengths is weaker 
than that at the P-branch wavelengths for a given 
3u by the same Boltzmann factor that the lower 
level population is greater, resulting in approxi- 
mately the same number of upward transitions in 
the two lines. The temperature gradient in the 
dust shell model caused the color temperature of 
the continuum radiation field to be higher than 
the gas temperature where the lines are formed. 
As a result, heat must flow from the radiation field 
into the gas. This happens by net absorption of 
R-branch photons and net emission of P-branch 
photons, pumping the rotational populations. 

It would be very desirable to estimate the mag- 
nitude of the error that is made using the isother- 
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mal absorbing slab model. For the models shown 
here, R-branch lines are typically a factor ~2 
weaker when reemission is included, but P-branch 
lines may be dominated by emission. Unfortu- 
nately, the extent to which emission cancels ab- 
sorption depends on the probability that an emit- 
ted photon escapes, rather than being absorbed 
by dust, and the asymmetry of the source, which 
determines how the absorption along our line of 
sight compares to that along other lines of sight. 
An understanding of the geometry of the source 
being observed would be needed to make a realis- 
tic model including reemission. 

Several simplifications and approximations 
made in the model should be noted. First, the 
gas temperature was assumed to be equal to the 
dust temperature. This assumption may be valid 
through much of the modeled shell, even if the 
density is too low for collisions with dust grains 
to dominate the heating and cooling of the gas, 
since the radiation field temperature is close to 
the dust temperature through most of the shell. 
However, close to the exciting star ultraviolet ra- 
diation may heat the gas to a higher temperature 
than the dust, via the photoelectric effect. The 
optical depth of the dust in the model is high 
enough to hide this region, but in the case of less 
deeply embedded objects, where the gas near the 
exciting source can be observed, departure of the 
gas temperature from the dust temperature could 
be important. To test this effect, a less optically 
thick model was run that included photoelectric 
heating of the gas and allowed the gas and dust 
temperature to differ. In this model, the gas tem- 
perature was greater than the infrared radiation 
temperature, so energy should flow from the gas 
into the radiation. In fact, P-branch lines then 
showed greater absorption than R-branch lines, as 
is expected by this thermodynamic argument. 

Another assumption was that the rotational 
temperature of the molecules was equal to the gas 
temperature. This might be justified by the fact 
that the density in most of the shell is greater than 
the critical density for rotational thermalization, 
especially for symmetric molecules, like C2II2 and 
CH4, which have no allowed rotational transitions. 
However, the infrared pumping that led to differ- 
ing P and R-branch line proflles may modify the 
populations. The importance of infrared pumping 
can be estimated by comparing the rate of coUi- 



sional transitions between rotational levels to the 
rate of infrared transitions between the ground 
and excited vibrational levels. The collisional 
transition rate is given by the product of the colli- 
sional rate coefRcient, which is typically ~ 10~^°, 
and the gas density, giving a rate ~ 2 x 10~^s~^. 
The radiative rate is given by the product of the 
vibrational A coefRcient and the fractional popu- 
lation of the excited vibrational state. For C2H2, 
with an A coefficient of 6 s~^ for the sum of the 
P and R-branch lines to a given rotational state, 
infrared pumping should dominate for vibrational 
temperatures > 130 K, or throughout much of the 
model shell. This calculation may overestimate 
the importance of infrared pumping since absorp- 
tion followed by emission of a vibration-rotation 
photon can at most change the rotational quan- 
tum number by 2, whereas large A J can occur in 
collisional transitions. But if the rotational pop- 
ulations are controlled by infrared pumping, the 
rotational temperature should equal the infrared 
radiation color temperature, and the difference be- 
tween P and R-branch lines should disappear. The 
resolution of this problem may be that the density 
of the gas around IRS 9 is enough greater than that 
assumed in the model so that collisions, rather 
than infrared pumping control the rotational level 
populations. A similar effect is the pumping of vi- 
brational levels by radiative and collisional transi- 
tions to other vibrational states. For example, the 
i/5 state of C2H2 could be populated by collisional 
or radiative transitions from the 1^4 state, which 
itself has several possible methods of excitation. 
To include the various neglected ways in which 
the excited vibration-rotation states can be pop- 
ulated would require the simultaneous calculation 
of many vibration-rotation states. Unfortunately, 
that is beyond the capability of the model used 
here. 

6. Conclusions 

The flrst conclusion of this work is that in- 
terpretation of vibration-rotation spectra based 
on the isothermal absorbing slab model can be 
very misleading. In general, column densities de- 
rived with this model can be expected to be un- 
derestimates, due to the neglect of reemission by 
molecules following absorption, but the magnitude 
of the error depends on the source geometry and 
the lines observed. 
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The second conclusion is more positive. Al- 
though collisions are not normally sufficiently fre- 
quent to maintain LTE at the gas kinetic tem- 
perature, the radiation field temperature may be 
similar to the gas temperature, and so molecular 
level populations may nevertheless be close to LTE 
populations. This was the case for the model con- 
sidered, and should generally be the case if the 
dust optical depth is large. As a result, non-LTE 
effects may not be large. Effects of radiative trans- 
fer through gas and dust with a temperature gra- 
dient are likely to be more important. Fortunately, 
these effects arc easier to include in models used 
to interpret observations, and they should be in- 
cluded. 

I thank the anonymous referee for several in- 
teresting comments, in particular pointing out the 
possibility of infrared pumping of the rotational 
populations. 
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